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Abstract 

We study the properties of the Arnold cap map on a torus with a several periodic sections using the Ulam method. 
This approach generates a Markov chain with the Ulam matrix approximant. We study numerically the spectrum and 
eigenstates of this matrix showing their relation with the Fokker-Plank relaxation and the Kolmogorov- Sinai entropy. 
We show that, in the frame of the Ulam method, the time reversal property of the map is preserved only on a short 
Ulam time which grows only logarithmically with the matrix size. Parallels with the evolution in a regime of quantum 
chaos are also discussed. 



1. Introduction 

The Arnold cat map flj] is the cornerstone model of 
classical dynamical chaos 12-0] • This symplectic map 
belongs to the class of Anosov systems, it has the pos- 
itive Kolmogorov- Sinai entropy h ^ 0.96 and is fully 
chaotic 1 4]. The map has the form 

^ = /? -h X (mod L) , x = x + p{mo&l). (1) 

Here the first equation can be seen as a kick which 
changes the momentum /? of a particle on a torus while 
the second one corresponds to a free phase rotation in 
the interval -0.5 < x < 0.5; bars mark the new values 
of canonical variables {x,p). The map dynamics takes 
place on a torus of integer length L in the p direction 
with -LI 2 < p < L/2. The usual case of the Arnold cap 
map corresponds to L = 1 but it is also possible to study 
the chaotic properties of the map on a torus of longer in- 
teger size L > 1 as it has been discussed e.g. in |5]. For 
L » 1 the spreading in p is characterized by a diff'usive 
process described by the Fokker-Planck equation: 



dw(p, t)/dt = D/2 d^w(p, t)/d^p. 



(2) 



where the diffusion coeflficient Z) ^< x^ >= 1/ 12 w(p,t) 
is a probability distribution over momentum and t be- 
ing integer time measured in number of iterations. As a 
result for times t » L^/D the distribution converges 
to the ergodic equilibrium with a homogeneous den- 
sity in the plane (x,p). The exponential convergence 



to the equilibrium state is determined by the second 
eigenvalue A2 of evolution Q on one map iteration with 
\A\ = Qxp(-To) < 1 and the convergence rate 



Yd = 271^ D I U- ^ 1M49/L^ 



(3) 
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the fist eigenvalue is Ai = 1. 

The dynamical equations ([T]) are reversible in time, 
e.g. at the middle of free rotation, but, due to chaos and 
exponential instability of motion, small round-off' errors 
break time reversal leading to an irreversible relaxation 
to the ergodic equilibrium L5J . 

In this work we investigate the transition from dy- 
namical behavior to statistical description using the 
Ulam method proposed in 1960 |6]. According to this 
method the whole phase space is covered by equidistant 
lattice (N = Np X Nx in our case). Then the transition 
probabilities from cell to cell are determined by propa- 
gating a large number of trajectories Ntr from one initial 
cell j to all other cells / after one iteration of the map 
(we used here A^^^ = 10^). In this way we generate the 
Markov chain ||7|] with a transition matrix 5"/^ = Nij/Ntr, 
where Ntj is the number of trajectories arrived from 
cell j to cell /. By construction we have 2/Ii ^ij = 1 
and thus the matrix S belongs to the class of Perron- 
Frobenius operators 111 0, 1?] . It is proven that for hy- 
perbolic maps in one and higher dimensions the Ulam 
method converges to the spectrum of continuous sys- 
tem [9-1 1]. At the same time it is known that in certain 
cases the Ulam method gives significant modifications 
of the spectrum compared to the case of the continuous 
Perron-Frobenius operators flu]. Indeed, for Hamil- 
tonian maps with divided phase space the spectrum is 
completely modified (see discussions in |ll2lll3|]) due to 
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Figure 1: Complex spectrum A of the Ulam matrix approximant S 
for the Arnold cat map. Top three rows: right column has Nx = 47 
and L = 3, 4, 8 (from top to bottom); left column has N^ = 43 and 
L = 3,4,8 (from top to bottom). Bottom row is for L = 4 with 
A^;, = 50(left) and A^^ = 51 (right). The total matrix size is A^ = A^^A^;, 
with Np = LNx and all eigenvalues are shown for each panel. Unit 
circle is shown in red. 



penetration of trajectories inside stability islands. From 
a physical view point the discretization corresponds to 
an effective noise in canonical variables which ampli- 
tude is equal to the cell size. Since an arbitrary small 
noise gives propagation of trajectories inside stability is- 
lands 1 4] the spectrum of the Ulam matrix approximant 
of size N in such a case differs from its continuous limit. 
A generalization of the Ulam method, based on one er- 
godic trajectory, allows to obtain a con verg ent spectrum 
for dynamics on a chaotic component \Vy\ . 

The majority of numerical studies with the Ulam 
method has been done for one-dimensional maps (see 
e.g. lll4ll5lll7l] ) but recently the studies were extended 
to the two-dimensional maps (see e.g. | 12, 13, 16, 181). 
In a certain respect the interest to such studies was gen- 
erated by similarities between properties of the Ulam 
matrix approximant for dynamical maps, which can be 
viewed as the Ulam networks, and the Google matrix of 



the World Wide Web as it is discussed in yjl [iTI] . For 
2D dissipative maps it was found that th e sp ectrum is 
characterized by the fractal Weyl law lll2lll8n . 

In a difference from the previous studies of the Ulam 
method in 2D maps here we choose the Arnold cat map 
on a torus of size L since it is fully chaotic, it has well 
defined diffusive relaxation to the ergodic state at large 
L, and it is time reversible. Thus the aim of this work 
is to understand the interplay of all these features in the 
frame of the Ulam method and the finite size Markov 
chain with the Ulam matrix approximant S generated 
by this method. 

The paper is composed as follows: in Section 2 we 
describe the properties of spectrum and eigenstates of 
the matrix S , the features of time reversal are analyzed 
in Section 3 and discussion of the results is presented in 
Section 4. 



2. Spectrum and eigenstates of the Ulam matrix ap- 
proximant 

The complex eigenvalues At and right eigenvectors 
ij/i of the Ulam matrix approximant S satisfy the equa- 
tion Sij/i = Aiij/i and are determined numerically by 
direct dioganalization. In agreement with the Perron- 
Frobenius theorem 1 8] the maximal eigenvalue is /li = 1 
with the corresponding eigenstate being real, nonneg- 
ative and homogeneously distributed over the whole 
phase space. 

The global distributions of eigenvalues Xt in the com- 
plex plane are shown in Fig. [T] for even and odd number 
of cells. Usually we keep A^^ = LNx to have exactly the 
same amount of cells in each of L sections of the con- 
tinuous map. The results show that 2D distributions are 
different for even and prime values of A^;^ (see Fig. [T]). 
For the even case /I- values are homogeneous inside a 
circle of a certain radius. For the odd case the distribu- 
tion has a form of a ring without eigenvalues at \A\ ^ 
(or with a smaller density at zero). The arithmetic prop- 
erties of the number of cells Nx and A^^ = LNx play a 
visible role. Thus for A^;^ = 47 we have a formation of 
star with 16 star rays while for A^;^ = 43 there are 44 
rays which are much less visible (for A^;^ = 37 we obtain 
a similar type of distribution with 38 rays). We obtain a 
similar type of ring spectrum also for A^;^ = 5 1 (bottom 
right panel in Fig.[T]). In the case when both A^;^ and A^^ 
are primes, e.g. Nx = 47, A^^ = 191, (and hence we 
have only approximate relation A^^ ^ LNx) the visibility 
of rays also decreases (data not shown). 

We expect that in the limit of large A^;^ and A^^ = LNx 
with fixed L the distribution will converge to a limiting 
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Figure 3: (color online) Spectral gap A = 1 - I/I2I as a function of L. 
Solid (black) and dashed (red) curves represent the cases of N^ = 23 
and Nx = 47 respectively; Np = LNx- Dotted (blue) line shows the 
theoretical dependence A oc L~^. 



Figure 2: (color online) Dependence of density of states dWjdy on 
the decay rate 7 for the Ulam matrix approximant of the cat map 
with L = 4. The number of cells Nx and Np are even-even on 
top panel, prime-even on middle panel, and prime-prime on bottom 
panel. The corresponding values are shown on legends with the nota- 
tion Nx 'xNp. The densities of states are normalized by the condition 

^^^ dWldydy = 1. 



one in agreement with the spirit of mathematical the- 
orems about the convergence of Ulam matrix approxi- 
mant for fully chaotic maps |9-11]. A confirmation of 
this is seen in Fig. [2] where the density distributions of 
eigenvalues dW/dy are shown as a function of the re- 
laxation rate y = -21n|/l|. Indeed, the density is es- 
sentially size independent showing two distinct distri- 
butions for even and odd values of Nx. We suppose that 
this diff'erence between two cases can be related to the 
eff'ect of discretization on the continuous map symme- 
try X -^ -X. The third type of size independent dis- 
tribution appears in the case of prime values of Nx and 
Np ^ LNx (see Fig. [2]) but in this case the diff'erence 
should be attributed to the fact that this discretization 
does not preserve exactly L identical classical segments 
of the continuous map. 

The maximum of the distribution dW/dy is located 
approximately aty = 2 ^ 2h corresponding to the value 
of the Kolmogorov- Sinai entropy h. Thus these y val- 
ues describe the process of exponential divergence of 
nearby trajectories and are related to the exponential 
correlations decay generated by chaotic dynamics. In 
addition to these values 7 ~ 1 there is also the value 
of A2 = exp(-72/2) which is positive and is very close 
to the unit value Ai = 1. It corresponds to the second 
eigenvalue of the Fokker-Plank equation describing dif- 



fusive relaxation to the ergodic steady- state. Indeed, the 
dependence of the gap A = 1 - /I2, shown in Fig. [3l is 
well in agreement with the dependence ([3]): a formal fit 
for 3 < L gives A oc 1/L^ with the exponent yu = 1.97 
for A^;, = 23 and 1 .94 for A^;, = 47 being close to the the- 
oretical value [1 = 2. The fit at the fixed exponent// = 2 
gives the numerical constant a in the relation A = a/L^ 
being 1.56+0.05 for A^^, = 23 and 1.33+0.05 for A^^, = 47 
that is close to the theoretical value a = 27i^D = 1.64.... 
A small deviation can be attributed to the eff'ect of finite 
size discretization. 

Two examples of eigenstates 1/^2 and j/^27 of the matrix 
S are shown in Fig. [4] (the numbering is done in a de- 
creasing order of \A\). According to the Fokker-Planck 
equation ^ we expect to have two double degener- 
ate values of A2 with the corresponding running wave 
eigenstates i//2 oc exp(±ik27ip/L) with ^ = +1 or their 
liner combination. The numerically found eigenstate i//2 
is real up to a numerical level of precision of matrix di- 
agonalization, corresponding to a real A2 value. It shows 
a certain amplitude oscillations along p but its main fea- 
ture is the sign change along x, which is absent in the 
equation (O. A dependence of eigenstates on x variable 
remains strongly visible and for other eigenstates (see 
Fig.©. For example, the state i//27 has density concen- 
tration at points x = and x = +0.5 corresponding to 
zero force point and discontinuity point respectively. 

These results show that the Fokker-Planck equation 
gives only a fist approximation for the statistical de- 
scription of dynamics of the Arnold cat map. In contrast 
to that the Ulam matrix approximant gives much more 
detailed description. The understanding of all features 
of this statistical description requires further more de- 
tailed studies. We note that according to (O there should 
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Figure 4: (color online) Eigenstates of the matrix S shown on the 
whole phase space for L = 4, Nx = 41 and Np = LNx. Top panel 
shows the eigenfunction il/2, which is real and has A2 = 0.917; middle 
and bottom panels represent respectively real and imaginary part of 
the eigenstate (/^27 corresponding to the eigenvalue /I27 - 0.516 + 
/0.507. In all panels red corresponds to positive, and blue to negative 
values. The phase space is rotated on 90 degrees clockwise so that 
positive values of p are on the right hand part of the plot. 



be a series of eigenstates with Ak = expi-Tok^)- We 
well resolved the case of ^ = 0, ± 1 but we found difficult 
to define accurately the corresponding higher values. It 
is possible that they enter rapidly in the dense balk re- 
gion of ring with many eigenvalues and become mixed 
with them. Probably larger values of L should be stud- 
ied to resolve such eigenvalues in a better way. We note 
that such a series of eigenvalues has been seen for the 
Chirikov standard map at the critical value of chaos pa- 
rameter where the diff'usion rate is relatively small and 
thus these eigenvalues are better separated from the balk 
region [131 . 

3. Time reversal features and the Ulam time 

Even if the exact dynamics is time reversible it be- 
comes easily broken by small errors due to dynamical 
chaos and exponential instability of motion (see e.g. 
iSlllSl]). In the quantum case the evolution is described 
by the linear Schrodinger equation that together with the 
uncertainty principle leads to stability of time reversal in 
respect to small errors iH, H l20] • Let's now study the 
time reversal in the frame of the Ulam method where 
the evolution is described by a linear matrix transfor- 
mation. For that we start from an initial line in the 




Figure 5: (color online) Time evolution of (p^) for an initial distribu- 
tion with p = obtained with the Ulam matrix approximant of the 
Arnold cat map at L = 8, A^p = LN^ and N^ = II (top-left panel), 51 
(top-right panel), 251 (bottom-left panel) and 751 (bottom-right pan- 
els). The time inversion is done at tr = 2,6, 10, 14. Yellow dashed 
line shows the theoretical prediction with (p^) = Dt and D = 1/12. 
The arrows in top-left panel shows the values of time a, b, c, d plotted 
inFig.[6l 



phase space at /? = with a homogeneous density in 
X, with w{pj = 0) = \INx and zero otherwise. Then 
we follow the evolution given by the matrix multiplica- 
tion w(/?, xj+l) = Sw(;?, X, t) where time t is measured 
in the number of map iterations t. The growth of the 
second moment {p^(t)} as a function of time is shown in 
Fig.O The second moment grows diff'usely with time in 
agreement with the Fokker-Planck equation ^. After tr 
iterations we perform time reversal by inverting all mo- 
menta p ^ -p. We see that after tr the second moment 
{p^(t)} starts to decrease during a certain time interval 
tu, where its value becomes minimal, and after that the 
diff'usion restarts again. This time tu is the Ulam time 
scale during which we have anti-diff'usive process which 
also describes relaxation in a vicinity of big fluctuations 

The spreading of probability in the whole phase space 
is displayed in Fig. [6] at difl'erent moments of time. 
These data show that time reversal gives a tempo- 
rary shrinking of the distribution followed by difl'usive 
spreading continued. Thus the time reversal is preserved 
only on the Ulam time scale tu- 

The strong efl'ects of exponential instability on time 
reversal breaking are also well seen in Fig. [71 Here, the 
initial image of the Arnold cat cannot be recovered even 
if the time reversal is done after only tr = 4 map iter- 
ations for the discretization level with A^;^ = 51. Only 
for a much finer discretization level with A^;^ = 751 the 
initial image is approximately recovered for tr = 4 but 
it degradates rapidly already at time reversal performed 
after tr = 5. This illustares the exponentially rapid 




Figure 6: (color online) Time evolution of a p = state with the Ulam 
matrix of the Arnold cat map at L = 8, A^j^ = 51 and Np = LN^ = 408. 
Panels show distribution in the whole phase space at the moments of 
time a, b, c, d marked by arrows in Fig.|5] (a) ait = 6 (top-left panel); 
(b)a.tt= 10 with p-inversion in t = 6 (top-right panel); (c): t = 12 
with p-inversion int = 6 (bottom-left panel); (d) Sitt = 12 without p- 
inversion (bottom-right panel). Probability density is shown by color 
with blue for zero density and red for maximal density on a given 
panel. The values of {p^} are marked in green arrows on top-right 
panel of Fig. \5\ 




Figure 7: (color online) Evolution of the Arnold cat image generated 
by the Ulam matrix approximant S with Nx = 51 (left colum) and 
Nx = 151 (right column) at L = 8, A^p = LNx. From top to bottom: 
initial image t = 0; image at the moment of time reversal t = tr = 4; 
image at the moment of return t = S; image at the return moment 
t = 10 when the time reversal is made after t = tr = 5 iterations. Color 
is proportional to density with blue for zero and red for maximum; 
only the central initial section -0.5 < x < 0.5, -0.5 < x < 0.5 is 
shown. 
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Figure 8: (color online) Dependence of the Ulam time ty with anti- 
diffusion on the number of cells Nx of the Ulam matrix approximant 
at L = 8 and Np = LNx- An initial state at momentum p = and 
homogeneous in x is evolved up to time /^ = f^ = 30 when a time 
reversal operation is applied. The dashed line shows the dependence 
tu =lnNx/h 



breaking of time reversal after the Ulam time tu- 

The dependence of the Ulam time tu, characterized 
by the anti-diffusion during time reversal process seen in 
Fig. [51 on the discretization scale Nx is shown in Fig.[8l 
We see that the results are well described by the depen- 
dence 



tu = lnNJh = \lnheff\/2h, heff=l/Nl 



(4) 



where he// = l/NJ can be considered as an effec- 
tive Planck constant which gives the area of discretized 
cells. In this form the Ulam time scale tu is simi- 
lar to the Ehrenfest time tE 1221] which appears in the 
semiclassical limit of systems of quantum chaos. In 
both cases the mechanism is related to the exponen- 
tial growth of a wave packet of minimal size with time 
heff exp(ht) due to which the packet size becomes com- 
parable with the whole system size after time tE or in 
our case after time tu- 

In spite of this similarity we should note that in the 
quantum case the time reversal is preserved under rather 
generic conditions (see iSllSlIZO] andRefs. therein). In 
contrast to that in the frame of the Ulam method, which 
also describes the linear matrix evolution, the time re- 
versal is broken after the Ulam time tu . The main reason 
of this difference is related to the fact that the Ulam ma- 
trix approximant describes dynamics with eigenvalue 
modulus smaller than unity while the quantum dynam- 
ics is unitary. This result can be also understood from 
the view point of noise which have a size of discretized 
cells and which also breaks time reversal. 



4. Discussion 

In this work we studied the properties of the Ulam 
matrix approximant S , generated by the Ulam method, 
for the Arnold cat map on a torus of a few integer sec- 
tions L. We show that the spectrum eigenvalues of 
S converges to a limiting distribution in the limit of 
small cell discretization and large matrix size N. The 
main part of this spectrum have relaxation rates y with 
approximate values of the Kolmogorov- Sinai entropy 
in this system. There are also eigenvalues with much 
smaller relaxation rate which is in a good agreement 
with the statistical description by the Fokker-Planck 
equation. 

The continuous model has the property of time re- 
versal but in the frame of the Ulam method the time 
reversibility is broken on the Ulam time scale tu which 
grows only logarithmically with the decrease of the cell 
size in the Ulam method. Such a dependence has cer- 
tain parallels with that one found for the Ehrenfest time 
scale in systems of quantum chaos. However, even 
if in both cases the evolution is described by the lin- 
ear matrix equations the quantum systems preserves the 
property of time reversal in presence of weak pertur- 
bations, while for the Ulam method the time reversal 
is broken after the Ulam time scale tu- Further stud- 
ies are required for a better understanding of relations 
between the spectrum of the Ulam matrix approximant, 
chaos, diffusion, correlations decay and other statisti- 
cal properties of dynamical chaos in the Arnold cat map 
and other chaos systems. Thus this simple model of 
Vladimir Arnold still keeps its scientific wonder. 

Acknowledgments: we thank Elena Chepelianskii 
for cat image drawing used here and in ll5lll9l]. 
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